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Abstract 

In this talk I first review at an elementary level a selection of cen- 
tral limit theorems, including some lesser known cases, for sums and 
maxima of uncorrelated and correlated random variables. I recall why 
several of them appear in physics. Next, I show that there is room for 
new versions of central limit theorems applicable to specific classes of 
problems. Finally, I argue that we have insufficient evidence that, as 
a consequence of such a theorem, g-Gaussians occupy a special place 
in statistical physics. 
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Central limit theorems play an important role in physics, and in particular 
in statistical physics. The reason is that this discipline deals almost always 
with a very large number N of variables, so that the limit N — > oo required 
in the mathematical limit theorems comes very close to being realized in 
physical reality. Before looking at some hard questions, let us make an 
inventory of a few things we know. 

1 Sums of random variables 

Gaussians and why they occur in real life 

Let p(x) be an arbitrary probability distribution of zero mean. Draw 
from it independently N variables x\, x 2 , ■ ■ ■ , and then ask what is the 
probability Pn(Y) that the scaled sum (x\ + x 2 + ■ ■ ■ + xn)/N 1 I 2 take the 
value Y. The answer, as we explain to our students, is obtained by doing 
the convolution Pn(Y) = p(xi) * p(x 2 ) * • • -p{xn)- After some elementary 
rewriting one gets 

where the dots stand for an infinite series of terms that depend on all mo- 
ments ofp(x) higher than the second one, (x 3 ), (x A ), .... In the limit N — > 00 
the miracle occurs: the dependence on these moments disappears from ([!]) 
and we find the Gaussian P 00 (y) = (2vr(x 2 )) 1 / 2 exp {-\Y 2 /(x 2 )) . 

The important point is that even if you didn't know beforehand about 
its existence, this Gaussian results automatically from any initially given 
p(x) - for example a binary distribution with equal probability for x = ±1. 
This is the Central Limit Theorem (CLT); it says that the Gaussian is an 
attractor pQ under addition of independent identically distributed random 
variables. An adapted version of the Central Limit Theorem remains true 
for sufficiently weakly correlated variables. 

This theorem of probability theory is, first of all, a mathematical truth. 
In order to see why it is relevant to real life, we have to examine the equations 
of physics. It appears that these couple their variables, in most cases, only 
over short distances and times, so that the variables are effectively indepen- 
dent. This is the principal reason for the ubiquitous occurrence of Gaussians 
in physics (Brownian motion) and beyond (coin tossing). Inversely, the pro- 
cedure of fitting a statistical curve by a Gaussian may be considered to have 
a theoretical basis if the quantity represented can be argued to arise from a 
large number of independent contributions, even if these cannot be explicitly 
identified. 
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Figure 1: The one-sided Levy distribution Li ^Y) is one possible outcome of a 
central limit theorem. 

Levy distributions 

Symmetric Levy distributions. Obviously the calculation leading to ([I]) 
requires that the variance (x 2 ) be finite. What if it isn't? That happens, in 
particular, when for x — > ±oo the distribution behaves as p(x) ~ c±|x| _1_a 
for some a £ (0, 2). 

Well, then there is a different central limit theorem. The attractor is a 
Levy distribution L a ^(Y), where Y is again the sum of the x n scaled with 
an appropriate power of N and where (3 £ (0, 1) depends on the asymmetry 
between the amplitudes c+ and c_. A description of the L a ^ is given, e.g., 
by Hughes ([2], see §4.2-4.3). In the symmetric case c + = c_ we have (3 = 0. 
Then L 1)0 (F) is the Lorentz-Cauchy distribution P^iY) = + Y 2 )] and 

L 2 fl(Y) is the Gaussian discussed above. 

Asymmetric Levy distributions. If c_ = 0, that is, if p(x) has a slow 
power law decay only for large positive x, then we have (3 = 1 and the limit 
distribution is the one-sided Levy distribution. In the special case a = |, 
shown in Fig. [I], we have the Smirnov distribution L\ v It has the explicit 
analytic expression 



All these Levy distributions are attractors under addition of random vari- 
ables, just like the Gaussian, and each has its own basin of attraction. 




(2) 



and decays for large Y as ~ Y 3 / 2 . 



3 




Figure 2: A sequence of probability distributions developing a "fat tail" propor- 
tional to ~ r~ 3//2 when k — > oo. 



Addition of nonidentical variables 

Mathematicians tell us that there do not exist any other attractors, 
at least not for sums of independent identically distributed (i.i.d.) vari- 
ables. However, suppose you add independent but non-identical variables. 
If they're not too non-identical, you still get Gauss and Levy distributions. 
The precise premises (the "Lindeberg condition" ) under which the sum of a 
large number of nonidentical variables is Gaussian distributed, may be found 
in a recent review for physicists by Clusel and Bertin [3] . 

Now consider a case in which the Lindeberg condition does not hold. 
Suppose a distribution poo(r), defined for r > 0, decays as r" 3//2 in the 
large-r limit. Let it be approximated, as shown in Fig. [21 by a sequence of 
truncated distributions Pi(t) , p 2 (r) , . . . ,p k (r), . . . which is such that in the 
limit of large k the distribution Pk{ T ) has its cutoff at r ~ k 2 . Then how will 
the sum ti = T\ + r 2 + . . . + tl be distributed? 

The answer is that it will be a bell-shaped distribution which is neither 
Gaussian nor asymmetric Levy, but something in between. It's given by a 
complicated integral that I will not show here. It is again an attractor: it 
does not depend on the shape of p{r) and the Pk(r) for finite r, but only 
on the asymptotic large-r behavior of these functions as well as on how the 
cutoff progresses for asymptotically large k. 

An example: support of ID simple random walk. An example of just 
these distributions occurs in the following not totally unrealistic situation, 
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Figure 3: A random walker on the positive half-line pays its first visit to site L at 
time ti- The sites (squares) already visited are colored black. The time interval 
between the first visits to L and to L + 1 is equal to T£+i. During this time 
interval the walker makes an excursion in the direction of the origin, as indicated 
by the dotted trajectory. The probability distributions pl(t~l) are independent 
but non-identical. 

depicted in Fig. El and which was studied by Hilhorst and Gomes [4]. A ran- 
dom walker on a one-dimensional lattice with reflecting boundary conditions 
in the origin visits site L for the first time at time t^. We can then write 
t L = Ti + t 2 + . . . + t Li where is the time difference between the first visit 
to the (k — l)th site and the kth site. As k increases, the tend to increase 
because of longer and longer excursions inside the region already visited. The 
Tfc are independent variables of the type described above. For L — > oo the 
probability distribution of ti tends to an asymmetric bell-shaped function of 
the scaling variable t^/L 2 which is neither Gaussian nor Levy. It is given by 
an integral that we will not present here. It is universal in the sense that it 
depends only on the asymptotic large r behavior of the functions involved. 

Sums having a random number of terms 

The game of summing variables still has other variations. We may, for 
example, sum N i.i.d. where A" itself is a random positive integer. Let 
A^ have a distribution tt^{v) where v is a continuous parameter such that 
(N) = v. Then for v — > oo one easily derives new variants of the Central 
Limit Theorem. 
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2 Maxima of random variables 

Gumbel distributions 



Let us again start from N independent identically distributed variables, 
but now ask a new question. Let there be a given a probability law p(x) which 
for large x decays faster than any power law (it might be a Gaussian). And 
suppose we draw N independent random values X\, . . . , x^ from this law. We 
will set Y = max!<j<7v %i- Then what is the probability distribution P^(Y) 
of Y in the limit Y — > oo ? The expression is easily written down as an 
integral, 



The calculation is a little harder to do than for the case of a sum. Let us 
subject the variable Y to an appropriate (and generally iV-dependent) shift 
and scaling and again call the result Y. Then one obtains 



which is the Gumbel distribution. 

The asymptotic decay of p(x) was supposed here faster than any power 
law. If it is as a power law, a different distribution appears, called Frechet; 
and if p(x) is strictly zero beyond some cutoff third distribution 

appears, called Weibull. Again, mathematicians tell us that for this new 
question these three cases exhaust all possibilities. 

In Ref. [3] an interesting connection is established between distributions 
of sums and of maxima. 

The Gumbel-k distribution. The Gumbel distribution (jlj) is depicted in 
Figs.Hl where it is called "Gumbel-1". This is because we may generalize the 
question and ask not how the largest one of the Xi, but how the kth largest 
one of them is distributed? The answer is that it is a Gumbel distribution 
of index k. Its analytic form is known and contains k as a parameter. For 
k — > oo it tends to a parabola, that is, to a Gaussian. 

All these distributions are attractors under the maximum operation. Even 
if you did not know them in advance, you would be led to them starting from 
an arbitrary given distribution p(x) within its basin of attraction. 

Bertin and Clusel (5j [6] show that the definition of the Gumbel-/c dis- 
tribution may be extended to real k. These authors also show how Gumbel 
distributions of arbitrary index k may be obtained as sums of correlated vari- 
ables. Their review article [3] is particularly interesting. 





Poo(Y)=e- Y e 



(4) 



The BHP distribution 
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Figure 4: The Gumbel-A: distribution for various values of k. 

In 1998 Bramwell, Holdsworth, and Pinton (BHP) [7J adopted a semi- 
empirical approach to the discovery of new universal distributions. These 
authors noticed that, within error bars, exactly the same probability distribu- 
tion is observed for (i) the experimentally measured power spectrum fluctu- 
ations of 3D turbulence; and (ii) the Monte Carlo simulated magnetization 
of a 2D XY model on an L x L lattice at temperature T < T c , in spin wave 
approximation. 

For the XY model Bramwell et al. [HI IS] were later able to calculate this 
distribution. It is given by a complicated integral that I will not reproduce 
here and is called since the "BHP distribution." Fig. shows it together with 
the Gumbel-1 distribution [TO] . 

Numerical simulations. How universal exactly is the BHP distribution? 
Bramwell et al. [8] were led to hypothesize that the BHP occurs whenever you 
look for the maximum of, not independent, but correlated variables. To test 
this hypothesis these authors generated a random vector X — (^1 i • • • 5 % N ) of 
N elements distributed independently according to an exponential, and acted 
on it with a fixed random matrix M such as to obtain y = Mx. By varying 
x for a single fixed M they obtained the distribution of Y = max!<j<jv y% 
and concluded that indeed it was BHP. 

However, Watkins et al. [TTJ showed one year later by an analytic calcu- 
lation that what appears to be a BHP distribution in reality crosses over to 
a Gumbel-1 law when N is increased. In this case, therefore, the correlation 
is irrelevant and the attractor distribution is as for independent variables. 

Watkins et al. conclude that "even though subsequent results may show 
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Figure 5: The Gumbel-1 and the BHP distribution. 



that the BHP curve can result from strong correlation, it need not. n This 
example illustrates the danger of trying to attribute an analytic expression 
to numerically obtained data. 

In later work Clusel and Bertin [3] present heuristic arguments tending 
to explain why distributions closely resembling the BHP distribution occur 
so often in physics. 

Wider occurrence of Gumbel and BHP? 

The Gumbel and BHP distribution have been advanced to fit curves in 
situations where their occurrence is not a priori expected. Two examples 
from the literature that appeared this month illustrate this. Palassini [12] 
performs Monte Carlo simulations that yield the ground state energy of the 
Sherrington-Kirkpatrick model; this author fits his data by a Gumbel-6 dis- 
tribution ([12], Fig. 4b). 

Gongalves and Pinto [32] consider the distribution of the cp daily return 
of two stock exchange indices (DJIA30 and S&P100) over a 21 year period. 
They find that the cubic root of the square of this distribution is extremely 
well fitted by the BHP curve ([13], Figs. 1 and 2). 

In both cases the authors are right to point out the quality of the fit. But 
these examples also show that having a very good fit doesn't mean you have 
a theoretical explanation. 



Periodic random walk 
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Figure 6: A random walk on a finite interval with periodic boundary conditions. 
The random variable x represents the maximum deviation of the walk from its 
interval average. The probability distribution P{x) was calculated by Majumdar 
and Comtet [14] . 

3 Correlated variables 

In addition to the example discussed above, we will provide here two further 
examples of how the maximum of a set of correlated random variables may 
be distributed. These examples will illustrate the diversity of the results that 
emerge. 

Airy distribution. Fig. [6] shows the trajectory of a one- dimensional ran- 
dom walker in a given time interval, subject to the condition that the starting 
point and end point coincide. The walker's positions on two different times 
are clearly correlated. Let x denote the maximum deviation (in absolute 
value) of the trajectory from its interval average. 

Majumdar and Comtet jH] were able to show that this maximum dis- 
tance is described by the Airy distribution (distinct from the well-known 
Airy function), which is a weighted sum of hypergeometric functions that I 
will not reproduce here. It is again universal: Schehr and Majumdar [T5] 
showed in analytic work, supported by numerical simulations, that this same 
distribution appears for a wide class of walks with short range steps. It turns 
out [H], however, that the distribution changes if the periodic boundary con- 
dition in time is replaced by free boundaries. This therefore puts a limit on 
the universality class [I]. 

Magnetization distribution of Ising 2D at criticality. We consider a finite 
L x L two-dimensional Ising model with a set of short-range interaction con- 
stants { Jfc}. Its magnetization (per spin) will be denoted M = N' 1 J2f=i s «> 
where N = L 2 and the Sj are the individual spins. We ask what the distribu- 
tion Pl{M) is exactly at the critical temperature T = T c . This distribution 
can be determined, in principle at least, by a renormalization calculation 
which in its final stage gives 



L-sP L (M) = V(m,{L- y * Ui }) 



(5) 
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Figure 7: Qualitative behavior of the distribution Pi(M) of the magnetization of 
the 2D Ising model on a periodic L x L lattice. The sharp peaks for T > T c and 
T <T C are Gaussians. For T = T c and under suitable scaling Pl(M) tends in the 
limit T —* oo to a double-peaked universal distribution V(m); see Eq. ([5]) and the 
accompanying text. 

where m = L^M and where {y^} is a set of positive fixed-point indices with 
corresponding scaling fields {ue} (i.e. the ue are nonlinear combinations of 
the Jk). In the limit L — > oo the dependence on these scaling fields disappears 
and we have, in obvious notation, that lini^oo L~sP L (M) = Vim). 

In Fig. [7] the distribution P L [M) is depicted qualitatively for L ^> 1 (it 
has two peaks!), together with the Gaussians that prevail when T ^ T c . The 
reason for M not being Gaussian distributed exactly at the critical point is 
that for T = T c the spin pair correlation does not have an exponential but 
rather a slow power law decay with distance: the spins are strongly correlated 
random variables. 

The similarity between Eq. (j5J) and Eq. (OQ) is not fortuitous: the coarse- 
graining of the magnetization which is implicit in renormalization, amounts 
effectively to an addition of spin variables; and the set of irrelevant scaling 
fields {ue} plays the same role as the set of higher moments {(x n ) \ n > 3} 
in Eq. (pQ). 

Eq. (JS]) says that V(m) is an attractor under the renormalization group 
flow; it is reached no matter what set of coupling constants { Jk} was given at 
the outset. Here, too, there are limits on the basin of attraction: the shape 
of V(m) depends, in particular, on the boundary conditions (periodic, free, 
or otherwise [16J). 

The conclusion from everything above is that attractor distributions come 
in all shapes and colors, and that it makes sense to try and discover new ones. 
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4 g-Gaussians 



A q- Gaussian G q (x) is the power of a Lorentzian, 



G q (x) 



est 



est 



(6) 



[1 + ax 2 } 



[1 + (q- l)x 2 }^ 



where in the second equality we have set p = l/(g — 1) and scaled x such 
that a = q — 1. Examples of g-Gaussians are shown in Fig.[BJ For q = 2 
the g-Gaussian is a Lorentzian; in the limit q — > 1 it reduces to the ordinary 
Gaussian; for q < 1 it is a function with compact support, defined only for 
— x m < x < x m where x m = l/i/l — q. For q = it is an arc of a parabola 
and for q — > — oo (with suitable rescaling of x) it tends to a rectangular block. 

Interest in g-Gaussians in connection with central limit theorems stems 
from the fact that they have many remarkable properties that generalize those 
of ordinary Gaussians. One may consider, for example, the multivariate q- 
Gaussian obtained by replacing x 2 in ((HD with v=l x^A^ y x v (with A a 
symmetric positive definite matrix). Upon integrating this g-Gaussian on 
m of its variables we find that the marginal (n — m)-variable distribution is 
g m -Gaussian with q m = 1 — 2(1 — g)/[2 + m(l — g)] (see Vignat and Plastino 
[17) ; this relation seems to have first appeared in Mendes and Tsallis [18]). 

A special case is the uniform probability distribution inside an n-dimensional 
sphere of radius R, 



where 6 denotes the Heaviside step function. This is actually a multivariate q- 
Gaussian with q = — oo. Integrating on m of its variables yields a g-Gaussian 
with q m = 1 — 2/m. We see that for large m both in the general and in the 
special approaches unity and hence these marginal distributions tend 

under iterated tracing to an ordinary Gaussian shape. 

Let us first see, now, how g-Gaussians may arise as solutions of certain 
partial differential equations in physics. 

Differential equations and g-Gaussians 

Thermal diffusion in a potential. The standard Fokker-Planck (FP) equa- 
tion describing a particle of coordinate x diffusing at a temperature T in a 
potential U(x) reads 







11 



Figure 8: Examples of q-Gaussians. Dotted curve: the ordinary Gaussian, q = 
1. Solid curves: the g-Gaussians for q = 2.2 and q = —2; the former has fat 
tails whereas the latter is confined to a compact support. All three curves are 
normalized to unity in the origin. 

Its stationary distribution Pfj(x) is the Boltzmann equilibrium in that poten- 
tial, Pfj(x) = est x exp[—/3U(x)], where (3 = 1/ksT. For the special choice 
of potential U'(x) = ax /{I + jx 2 ) the stationary distribution becomes the 
g-Gaussian 

P$(x)=cstx (l + 7 x 2 )- a/3/7 . (9) 

This distribution is an attractor under time evolution, the latter being defined 
by the FP equation (jSJ); a large class of reasonable initial distributions will 
tend to (jUJ) as t —> oo [19J. It should be noted, however, that by adjusting 
U(x) we may obtain any desired stationary distribution, and hence the q- 
Gaussian of Eq. (|5J) plays no exceptional role. 

The following observation is trivial but will be of interest later on in 
this talk. Let x(t) be the Brownian trajectory of the diffusing particle. 
Let x(0) be arbitrary and let x(t), for t > 0, be the stochastic solution 
of the Langevin equation associated [20] with the FP equation (jSj). Let 
£ n = x(nr) — x((n — l)r), where r is a finite time interval. Then Y/v = 
£i + ^2 + • • • + (,n (without any scaling) is a sum which for iV — >• oo has the 
distribution P^(Y). In particular, if U(x) is chosen such as to yield Qj, we 
have constructed a g-Gaussian distributed sum. 

Finite difference scheme [2TJ. Rodriguez et al. [22] recently studied the 
linear finite difference scheme 

TN,n + r N , n+ l = r N _ 1:U , (10) 

where N = 0, 1, 2, . . . and n = 0, 1, . . . , N. The quantity pn,u = ( n )r/v,n may 
be interpreted as the probability that a sum of N identical correlated binary 
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variables be equal to n. For specific boundary conditions, the authors were 
quite remarkably able to find a class of analytic solutions to Eq. (11 01) and 
observed that the N — > oo limit of the sum law pn,u is a g-Gaussian. 

To understand better what is happening here, let us set t = 1/N, x = 
1 — 2n/N, and P{x 1 t) = N( )r^ jn . When expanding Eq. (fIU|) in powers of 
iV -1 one discovers [23J that P(x,t) satisfies the Fokker-Planck equation 

^-&i*->W (id 

for t > and — 1 < x < 1. The "time" t runs in the direction of decreasing 
N. Hence Rodriguez et al. have solved a parabolic equation backward in 
time and determined, starting from the small- N behavior, what is actually 
an initial condition at N = oo. It is obvious that g-Gaussians are not singled 
out here: there exists a solution to Eq. (fTTT) for any other initial condition at 
t = 0, and concomitantly to Eq. (11 01) for any desired limit function Poo tn at 
N = oo. 

Therefore, in this and the preceding paragraph, the occurrence of q- 
Gaussians in connection with Fokker-Planck equations cannot be construed 
as an indication of a new central limit theorem. 



The porous medium equation. Let us consider a fluid flowing through a 
porous medium. Three equations of physics provide the basic input for the 
description of this flow, namely (i) the continuity equation for the fluid den- 
sity p(x,t); (ii) Darcy's law, which relates the fluid velocity v to its pressure 
p by v = —est x Vp; and (iii) the adiabatic equation of state of the ideal gas. 
Upon combining these one finds the porous medium equation 

^| = Ap 2 -«, q=l-C p /C v , (12) 

where C p /C v is the specific heat ratio. For q = 1 this equation reduces to 
the ordinary diffusion equation. 

Equation (1121) is nonlinear and its general solution, i.e., for an arbitrary 
initial condition u(x, t), cannot be found. It is however possible to find special 
classes of solutions. One special solution is obtained by looking for solutions 
that are (i) radially symmetric, i.e., dependent only on x = \x\; and (ii) scale 
as u(x,t) = t~ db F(xt~ b ). After scaling of x and t we obtain the similarity 
solution 



x 2 



1 + (1 ~ 1) t2l , 



(13) 



in which b = — g) + 2] and where also Cq is uniquely defined in terms 

of the parameters of the equation. Mathematicians (see e.g. [21]) have shown 
that initial distributions with compact support tend asymptotically towards 
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this similarity solution. The asymptotic behavior ({TBI) is conceivably robust, 
within a certain range, against various perturbations of the porous medium 
equation. It is not clear to me if and how this property can be connected to 
a central limit theorem. 

g-statistical mechanics 

Considerations from a (/-generalized statistical mechanics [251 [261 [27] have 
led Tsallis [28] to surmise that in the limit iV — > oo the sum of iV correlated 
random variables becomes, under appropriate conditions, g-Gaussian dis- 
tributed; that is, on this hypothesis g-Gaussians are attractors in a similar 
sense as ordinary Gaussians. Now, variables can be correlated in very many 
ways. To fully describe N correlated random variables you need the N vari- 
able distribution P N (xi, . . . , x^). Taking the limit iV — > oo requires knowing 
the set of functions 

P N ( Xl ,...,x N ), N = 1,2,3,... (14) 

In physical systems the P/v are determined by the laws of nature; the relative 
spatial and/or temporal coordinates of the variables, usually play an essen- 
tial role. The examples of the Ising model and of the Airy distribution show 
how widely the probability distributions of strongly correlated variables may 
vary. Hence, in the absence of any elements of knowledge about the physical 
system that they describe, statements of uniform validity about correlated 
variables cannot be expected to be very specific. 

g-Central Limit Theorem 

We now turn to a g-generalized central limit theorem (g-CLT) formulated 
by Umarov et al. [25]. It says, essentially, the following. Given an infinite set 
of random variables xi, X2, ■ ■ ■ , x n , . . . , let the first iV of them be correlated 
according to a certain condition Cjv(g), where N = 1,2,3,.... Then the 
partial sum Yn = J2n=i x m after appropriate scaling and in the limit iV — ► 
oo, is distributed according to a g-Gaussian. The theorem is restricted to 1 < 
q < 2. The conditions C^{q) are concisely referred to as "g-independence" 
in Ref. [2H] and for g = 1 reduce to the usual condition of random variables 
being independent. Closer inspection of the theorem prompts two questions. 

First, the conditions C^iq) are difficult to handle analytically. If a the- 
oretical model is defined by means of its Pn&i, ■ ■ ■ , %n) f° r N = 1, 2, 3, . . ., 
then one would have to check that these satisfy the C^{q). I am not aware of 
cases for which this has been possible. In the absence of examples it is hard 
to see why nature would generate exactly this type of correlations among its 
variables. 
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Secondly, the proof of the theorem makes use of "g- Fourier space," the g- 
Fourier transform (g-FT) having been defined in Ref. [29J as a generalization 
of the ordinary FT. The g-FT has the feature that when applied to a q- 
Gaussian it yields a g'-Gaussian with q' = (1 + g)/(3 — g), for 1 < g < 3. 
Now the g-FT is a nonlinear mapping which appears not to have an inverse 
[30J. It is therefore unclear at present how the statements of the theorem 
derived in g-Fourier space can be translated back in a unique way to "real" 
space. 

5 The search for g-Gaussians 

Mean-field models 

Independently of this g-CLT Thistleton et al. |31j (see also Ref. [32]) 
attempted to see a g-Gaussian arise in a numerical experiment. These au- 
thors defined a system of iV variables Xi, i = 1, 2, . . . , N, equivalent under 
permutation. Each variable is drawn from a uniform distribution on the in- 
terval (— f) i) but ^he Xi are correlated in such a way that (xjXk) = p(x\) 
for all j k, where p is a parameter in (0, 1) [33J. They considered the sum 
Yn — (xi + . . . + xn)/N and determined its distribution P(Y) in the limit 
N ^> 1. For p = jq the numerical results for P(Y) can be fitted very well by 
a g-Gaussian G q (Y) with g = — |, shown as the dotted curve in Fig.0 This 
system of correlated variables is sufficiently simple that Hilhorst and Schehr 
[31] were able to do the analytic calculation of the distribution. They found 
that Y is distributed according to 

P(y) = (2Z£) 4 exp (_ 2 SLlA [ erf -i(2y) ] 2 ) , (15) 

for — i < Y < i shown as the solid curve in Fig.0 The difference between the 
exact curve and the g-Gaussian approximation is of the order of the thickness 
of the lines. More importantly, the calculation of Ref. [31] shows that the 
distribution of the sum Y varies with the initially given one of the Xj. This 
initial distribution may be fine-tuned such as to lead for iV —>■ oo to almost 
any limit function P(Y) - in particular, to a g-Gaussian. The existence of 
g-Gaussian distributed sums was already pointed out below Eq. (jUJ) and is no 
surprise. However, there is, here no more than in the case of the FP equation, 
any indication that distinguishes g-Gaussians from other functions. 

The work discussed here concern a mean-field type model: there is full 
permutational symmetry between all variables. This will be different in the 
last two models that we will now take a look at. 

Logistic map and HMF model 
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Figure 9: Comparison of the g-Gaussian G q {Y) (dotted curve) guessed in Ref. [31] 
on the basis of numerical data and the exact distribution P(Y) (Eq. (I15h . solid 
curce) calculated in Ref. [33]. The curves are for p = ^ and the g-Gaussian has 
q = — §. The difference between the two curves is of the order of the thickness of 
the lines and just barely visible to the eye. 



Two well-known models of statistical physics have been evoked several 
times by participants [35j EE] at this meeting. The common feature is that 
in each of them the variable studied is obtained as an average along a deter- 
ministic trajectory. 

Logistic map. In their search for occurrences of g-Gaussians in nature, 
Tirnakli et al. [37] considered the logistic map 

xi = a — x\_ x , £ = 1,2,... (16) 

A motivation for this choice is the appearance [3S] of q- exponentials in the 
study of this map. Starting from a uniformly random initial condition x = Xq, 
Tirnakli et al. determined the probability distribution of the sum 

n a +N 

e=n 

of successive iterates, scaled with an appropriate power of N, in the limit 
N ^> 1. Their initial report of g-Gaussian behavior at the Feigenbaum crit- 
ical point (denned by a critical value a = a c ) was critized by Grassberger 
[39] . Inspired by a detailed study due to Robledo and Mayano [ID], who 
connect properties at a = a c to properties observed on approaching this crit- 
ical point, Tirnakli et al. [H] took a renewed look at the same question and 
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now see indications for a g-Gaussian distribution of Y near the critical point. 

Hamiltonian Mean Field Model. The Hamiltonian mean-field model (HMF), 
introduced in 1995 by Antoni and Ruffo |42j, describes L unit masses that 
move on a circle subject to a mean field potential. The Hamiltonian is, 
explicitly, 

w = Ef + ^ E [* - - > ( 18 ) 

i i,j 

where Pi and $i are the momentum and the polar angle, respectively, of the 
zth mass. The angles were originally considered to describe the state of 
classical XY spins, so that fhi = (cos 9{, sin 9j) is the magnetization of the 
ith spin. The HMF has a solvable equilibrium state. At a critical value 
U — U c = 0.75 of the total energy per particle a phase transition occurs 
from a high-temperature state with uniformly distributed particles to a low- 
temperature one with a spontaneous value of the "magnetization" (|M|), 
where M = L^ 1 ^2 i=1 fhi. When launched with certain nonequilibrium ini- 
tial conditions, the system, before relaxing to equilibrium, appears to enter 
a "quasi-stationary state" (QSS) whose lifetime diverges with N. It is im- 
possible to discuss here all the good work that has been done, and is still 
going on, to attempt to explain the properties of this state (see e.g. Chavanis 
H31 SI SSJ, Tsalhs et al. [46J, Antoniazzi et al. [47J, Chavanis et al. jig]). 
One specific type of numerical simulations, performed by different groups of 
authors, is relevant for this talk. These have been performed at the subcrit- 
ical energy U = 0.69 with initially all particles located at the same point 
(8i = for all i) and the momenta pi distributed randomly and uniformly in 
an interval [— p max , Pmax] ■ The QSS subsequent to these initial conditions has 
many features (such as non-Gaussian single-particle velocity distributions) 
that have been connected to g-statistical mechanics. Of fairly recent inter- 
est is the sum Yi of the single-particle momentum Pi(t) sampled at regularly 
spaced times t — It along its trajectory, 

n +N 

(19) 

£=n 

The distribution of Yi in the limit of large iV is again controversial [49, 50] • 
For the specific initial conditions cited above it seems to first approach a 
fat-tailed distribution, interpreted by some as a g-Gaussian, before it finally 
tends to an ordinary Gaussian. 

Comments. The analogy between (fTTj) and f[T9"|) is obvious. In both cases 
the sequence of iterates has long-ranged correlations in the "time" variable t 
and fills phase space in a lacunary way. It is therefore not very surprising that 
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Y and Yi should have non-Gaussian distributions. The q- Gaussian shape of 
these distributions, however, remains speculative. The examples of this talk 
have shown, on the contrary, that in the absence of specific arguments sums 
of correlated variables may have a wide variety of distributions. It seems 
unlikely that haphazard trials will hit exactly on the g-Gaussian. 



6 Conclusion 

Universal probability laws occur all around in physics and mathematics, and 
the quest for them is legitimate and interesting. What lessons can we draw 
from what precedes? 

• It is quite conceivable that new universal distributions may be dis- 
covered, either by asking new questions about independent variables; or by 
asking the traditional questions (sums, maxima,. . . ) about correlated vari- 
ables. 

• Variables may be correlated in an infinity of ways. In the end some 
real-world input is desirable, be it from physics, finance, or elsewhere. 

• Nothing can beat a central limit theorem. A good one, however, should 
give rise to analytic examples and/or simulation models that reproduce the 
theorem with high numerical precision. 

• In the absence of theoretical arguments, assigning analytic expressions 
to numerically obtained curves is a risky undertaking. 

Let me end by a quotation |51j: "Good theory thrives on reasoned dissent, 
and [our views] may change in the face of new evidence and further thought." 
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